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ABSTRACT 

The spectrum of the supernova relic neutrino (SRN) background from past stellar collapses includ¬ 
ing black hole formation (failed supernovae) is calculated. The redsliift dependence of the black hole 
formation rate is considered on the basis of the metallicity evolution of galaxies. Assuming the mass 
and metallicity ranges of failed supernova progenitors, their contribution to SRNs is quantitatively 
estimated for the first time. Using this model, the dependences of SRNs on the cosmic star formation 
rate density, shock revival time and equation of state are investigated. The shock revival time is intro¬ 
duced as a parameter that should depend on the still unknown explosion mechanism of core collapse 
supernovae. The dependence on equation of state is considered for failed supernovae, whose collapse 
dynamics and neutrino emission are certainly affected. It is found that the low-energy spectrum of 
SRNs is mainly determined by the cosmic star formation rate density. These low-energy events will 
be observed in the Super-Kamiokande experiment with gadolinium-loaded water. 

Subject headings: diffuse radiation — galaxies: evolution — neutrinos — supernovae: general 


1. INTRODUCTION 

Since the creation of the Universe, many generations of 
stars have been born and died. During the cosmic evolu¬ 
tion, stars eject synthesized elements by stellar winds or 
explosions such as supernovae, and the ejecta are mixed 
with the interstellar gas. Therefore, the mass fraction of 
elements heavier than carbon (metallicity), Z, increases 
gradually with the cosmic time. Meanwhile, many neu¬ 
trinos are emitted from core collapse supernova (CCSN) 
explosions of massive stars and accumulate to give a dif¬ 
fuse background radiation that is redshifted owing to cos¬ 
mic expansion. These neutrinos are called the supernova 
relic neutrino (SRN) background, or the diffuse super¬ 
nova neutrino background (DSNB) in some papers. 

Neutrinos emitted from a super nova have ac tually 
been detected for SN1987A (e.g . . iHirata et al.l 119871 : 
iBionta et al.l 119871 : lAlexevev et al.1119881) . In the obser¬ 
vation of SRNs, on the other hand, terrestrial neutrino 
detectors are affected by various backgrounds such as 
solar neutrinos, reactor neutrinos, atmospheric neutri¬ 
nos and contamination by cosmic muon events, radio 
activity events and so forth. However, some observa¬ 
tional upper bou nds fo r the flux of SRNs have been 
reported fe.g.. [Malek et al.112003 ). Roughly speaking, 
all species of neutrinos (z/ e , u e , u T , v T ) are 

equally emitted from a supernova with average energies 
of ~10 MeV. Nowadays, SRNs with v e of approximately 
20 MeV are expected to be observable in running exper¬ 
iments. The most stringent limits reported for v e flux 
were obtained in the Super-Kamiokande experiment as 
<0.1-1 cm -2 s" 1 MeV -1 f or neutr ino energ ies between 
17.3 MeV and 30.8 MeV ( Bavs et ah] 12012 ) and in the 
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Fig. 1.— 90% C.L. differential upper limits on v e flux of SRNs. 
The squares, circles and triangles are results for Super-Kamiokande 
rSK-I/II/III. IBays et al!2012fi. Sup er-Kamiokande with a neutron¬ 
taggi ng (SK-IV, Izhang et al .1120151) and KamLAND UGando et al.l 
2012). Dashed and dotted lines correspond to our theoretical mod¬ 
els with maximum and minimum values of SRN event rate, respec¬ 
tively (see also Table is. 

KamLAND experiment as <10-100 cm ~ 2 s " 1 MeV" 1 
between 8.3 MeV and 18.3 MeV ( Gando et al. 2012 ). 
Super-Kamiokande derived a new upper limit of <5- 
30 cm -2 s" 1 MeV" 1 for energies between 13.3 MeV 
and 17.3 MeV by performing a new analys is with a 
neutron-tagging technique (iZhang et al.l 120151) . In Fig¬ 
ure [Tj we show the upper limits for v & flux with the¬ 
oretical estimations presented later in this paper. For 
v e flux, the SNO experiment obtained an upper limit of 
70 cm" 2 s" 1 M eV" 1 f or ene rgies between 22.9 MeV and 
36.9 MeV (lAharmim et al.l I2006D . These observational 
upper limi ts are larger than various theoreti cal predic¬ 
tions fe.g.. lAndo fe Satoll2004 lBeacomll2010l and refer¬ 
ences therein). Nevertheless, the Super-Kamiokande up¬ 
per limit is reasonably close to the predictions; thus, it is 
expected that SRNs will be observed in the near future. 

Cosmic metallicity evol ution has been proven by obser¬ 
vations of galaxies (e.g., iMaiolino et al.l f2008. hereafter 
M08). Recently, the correlation between the metallic¬ 
ity and the star formation rate (SFR) of galaxies has 
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been studied for various ranges of the cosmic redshift pa¬ 
rameter z (e.g.. iMannucci et al.l 120101 : lYabe et all I2012L 
12014 iZahid et al.ll2014D . Since metals have strong effects 
on opacity and control the cooling rate and luminosity 
of star-forming clouds, SFR depends on the metallicity. 
On the other hand, SFR determines the amount of met¬ 
als ejected from stars and governs the metallicity evolu¬ 
tion. Among the newly born stars, those with mass of 
> 10M S undergo gravitational collapse to cause a CCSN 
or, if they fail to explode, form a black hole. Combining 
the SFR of individual galaxies with the galaxy distribu¬ 
tion, the cosmic star formation rate density (CSFRD) 
and total core-collapse rate can be determined at each 
redshift. Therefore, the flux of SRNs is closely related to 
the evolution of galaxies (|Totani et al.lll996l h 

Metallicity affects not only star formation but also stel¬ 
lar evolution. Since the stellar opacity and mass loss 
rate vary with the metallicity, the final fate of stars 
depends on the metallicity as well as the initial mass, 
M. Some progenitors may result in black-hole-forming 
failed supernovae depending on their mass and metallic¬ 
ity. They will contribute to the flux of SRNs because 
neutrinos are also emitted from them with higher lumi¬ 
nosit y and mean energy tha n those from ordinary CCSN 
(’e.g.- ISumivoshi et al.ll2006l ). The spectrum of SRNs that 
will be observed is the convolution of the neutrino spec¬ 
tra emitted from core collapses of progenitors with vari¬ 
ous masses and metallicities. Conversely, the observation 
of SRNs will provide useful information not only on su¬ 
pernova neutrinos themselves but also on cosmological 
evolutions of the star formation and metallicity. 

In this study, we construct a model spectrum of SRNs 
using the metallicity distribution function derived from 
models of galaxy evolution for the first time. Using 
this model, we investigate some uncertainties in SRNs, 
namely CSFRD, the explosion mechanism and the nu¬ 
clear equation of state (EOS). It is pointed out that most 
of the theoretical models for galaxy formation underpre¬ 
dict CSFRD compared with observations at z > 1 (e.g., 
iKobavashi et al.l 1 201. 'll hereafter K13). The explosion 
mechanism of a CCSN is still uncertain. Most numerical 
simulations reveal that the shock wave launched at the 
bounce of the inner core stalls on the way to the core sur¬ 
face and that some mechanism such as neutrino heating 
acts to revive the sho ck wave lead ing to the supernova 
explosion (e.g., iKotake et al.l 12012 ). The nuclear EOS, 
which is also still unknown, affects the collapse dynam¬ 
ics a nd neutrino emission, es pecially for failed supernovae 
(e.g., iSumivoshi et al.l[200611 . 

This paper is organized as follows. In § [2] we present 
the formulation of the SRN spectrum. Issues on neutrino 
oscillation are also given. In § [3] the models of galaxy 
evolution used in this study are described. We explain 
numerical models of the neutrino spectra emitted from 
various progenitors in § U Here we introduce the shock 
revival time as a parameter that reflects the unknown ex¬ 
plosion mechanism. In addition, we report a new result 
for a failed supernova with a different EOS. In § 0 we 
show the results for the spectrum of SRNs and the event 
rate for Super-Kamiokande over 1 year. We also investi¬ 
gate the uncertainties of the CSFRD, shock revival time 
and EOS for black hole formation. Finally, §H]is devoted 
to a conclusion and discussion. 


2. FORMULATION OF SUPERNOVA RELIC NEUTRINO 
BACKGROUND 

In this study, we construct a model spectrum of SRNs 
taking into account cosmic metallicity evolution. Since 
SRNs originate from various progenitors, their depen¬ 
dence of the neutrino emission affects the spectrum. Here 
we assume that the progenitors are characterized by their 
initial mass, M, and metallicity, Z. Then the flux of 
SRNs on the Earth is written as 


d F(E V ) 
d E v 


dz 


o F/q \ j H m (1 + z)^ + Ha 


Rcc{z) / *PzF(z,Z) 
Jo 




d E' 


( 1 ) 


with velocity of light c and cosmological constants Hq = 
70 km s -1 Mpc -1 , fl m = 0.3 and Ha = 0.7. The neu¬ 
trino energy on the Earth, E v , is related to that at the 
redshift 2 , E ' v , as E' v = (1 + z)E v . The total core¬ 
collapse rate, Rc c(z), is determined by CSFRD, which 
is stated in § 13.11 The neutrino number spectrum from 
the core collapse of a progenitor with mass M and metal¬ 
licity Z , dN{M 1 Z,E' v )/dE' l ,, includes the effect of neu¬ 
trino oscillation, which is stated below. Note that, in 
our model, not only ordinary supernovae but also black- 
hole-forming collapses without an explosion are consid¬ 
ered. Meanwhile, tp imf(AT) and -j/’zf {z,Z) are the ini¬ 
tial mass function and metallicity distribution function 
of progenitors, respectively, which are normalized as 
MF(M)dM : f 0 Z ™ x if> ZF (z,Z)dZ = 1. While 
the initial- mass dependence of progenitors was studied 
previously dNakazato et ahll2006HLunardini fe Tamborral 

2012; Nak azatoll2013H . the metallicity dependence is also 
considered in our study. Note that ip Z F (z, Z) is a function 
of the redshift z due to the cosmic metallicity evolution. 
We use the metallicity distribution function derived from 
models of galaxy evolution later in § 13.21 
As already mentioned, we should take into account 
neutrino oscillation in the evaluation of SRN spectra. 
Hereafter, with the Super-Kamiokande experiment in 
mind, we concentrate on v e . The number spectrum of 
v e in Equation m, d Np e (M, Z,E U )I d E v , is a mixture of 
spectra of neutrinos originally produced as v e , iy,, and 
v T and is written as 


dJVp.(M, Z,E„) 
d E v 

- d Nl(M,Z,E v ) t - dNl(M,Z,E„) 

— -Lp.p. , i ± ue 


d E v 

d Nl(M,Z,E„) 

i -l re 


d E„ 


= P 


d E v 

-diV? (M, Z,E V ) 


d E„ 


, _ d Nj> (. M,Z,E V ) , N 

+ (1~^) d E v ,(2) 


where P ae (a = e, //, r) is a conversion probability from 
v a to D e satisfying P ee + P^ + P Te = 1, and P (= P ee ) 
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is the survival probab ility of v e , passing through stellar 
envelopes and space (iDighe fc Smirnov! I2000f) . Mean¬ 
while, dlV? (M, Z, E v )/dE v is the original spectrum of 
D a . Note tliat the original spectra of and v T are al¬ 
most identical because and v T do not have charged- 
current reactions around the neutrino-emitting surface 
of a supernova core and, therefore, both of them can be 
denoted as dN® x (M, Z, E v )/dE v . These original spec- 
tra are taken from t he Supernova Neutrino Database 
dNakazato et al.ll2013lh 

The flavor eigenstates of antineutrinos, \n a ), are lin¬ 
ear combinations of the energy eigenstates in vacuum, 
\Di) (i = 1,2,3 with neutrino masses mi, m 2 ,m 3 , respec¬ 
tively): 

i^) = E^i p ‘)> ( 3 ) 

i 


inverted mass hierarchy, Equation 0 becomes 

d N*. d Nl 

dE v d E v ’ W 

where we neglect the small value of sin 2 9i 3 ~ 0.02. This 
means a complete transition (P = 0). Incidentally, we 
do not take into account the breaking of adiabati city 
due to the shock wave propagation (e.g., iGalais et al.l 
120101 ). which should be m inor for a time-integrated sig- 
nal dKawagoe et all I2010H . We also neglect neutrino- 
neutrino collect i ve eff ects (e.g.-IChakrab ortv et al.ll20Tlt 
iDasgnota et al.l l2012t) , whose c ontribution to SRNs is 
estim ated to be about 5-10% (|Lunardini fc Tamborral 

[M2). 

3. MODELS OF GALAXY EVOLUTION 


where the mixing matrix U can be expressed as 


Urvi — 


( C12C13 S12C13 

— S12C23 — Ci2S23Sl3e i<5 C12C23 — Sl2S23Sl3e I<5 
S12S23 ~ Ci2C23Si3e*' 5 —C12S23 — Sl2C23Sl3e Z<5 


si 3 e- a \ 
S 23 C 13 I 
C23C13 / 


( 4 ) 


with mixing angles 612 , 623 , 613 , CP phase S (we set S = 0 
in this study) and Sij = sin 0^, = cos 6 ^ . The 

three mixing angles have been measured in several ex¬ 
periments/observations: sin 2 0 i 2 ~ 0.31, sin 2 0 23 ~ 0.44, 
sin 2 0 i 3 ~ 0.02 (|Qlive et al.ll2014lh Using the matrix el- 
ements of U , the numb er spectrum of D e is denoted as 
(|Dighe fc Smirnovll2000D 


|tt 12 diV^j^ 
■ — |Uel| 


d E v 


d E v 


ITT | 2 dlW 2 

|c/e21 dP7 


I u, 12 


' e 3 | 


d N 03 
dE v 


2 n 2 n dlVy,^ .2/1 2 a dlV ^ 2 

= COS 012 COS 013——:;-b Sin 0 i2 COS 013- 


d E v 


d E v 


. 2 n dNp 3 
■sm 2 013-7T^: 
d E„ 


( 5 ) 


where dN^JdE^ is the number spectrum of z/;. 

In the supernova core, a good approximation is that 
the lowest neutrino energy eigenstate corresponds to v e 
and other states are mixtures of v M and v T . Furthermore, 
if the stellar density profiles are not sufficiently steep 
to break up the adiabaticity of propagating eigenstates, 
neutrinos emitted from the dense core pass through the 
stellar envelope with varying density along a single en¬ 
ergy eigenstate into vacuum space. Note that there is 
still uncertainty concerning the mass hierarchy: normal 
(mi < m 2 < m 3 ) or inverted (m 3 < mi < m 2 ). Thus, we 
obtain dN Dl /dE v ~ dN® e /dE v , dN^/dE v ~ dN® x /dE v , 
d-/Vp 3 /d E u ~ dNp x /dE v for the case of the normal mass 
hierarchy and dA % 3 /dE v ~ dN® e /dE v , dNp 1 /dE u ~ 
dN§ x /dE v , dNp 2 /dE l/ ~ dNp~/dE v for the case of the 
inverted mass hierarchy (|Dighe & Smirnov! 1200 0). Sub¬ 
stituting these relations into Equation (0, we obtain 


d N, e 
d E v 




( 6 ) 


and P = 0.68 for the normal mass hierarchy. For the 


To derive a realistic spectrum of SRNs, it is manda¬ 
tory to use a reliable core-collapse rate based on mod¬ 
els of galaxy evolution. Here we consider the stellar 
mass, M*, SFR, M*, and metallicity, 12+logio(0/H), 
of galaxies. The SFR of galaxies is responsible for 
the CSFRD, p*{z), and the metallicity of galaxies cor¬ 
responds to that of the progenitors in them. Fur¬ 
thermore, M*, M, and 12+logm(0/H) are related to 
each o ther (e.g..lBrin chmaniiet al.ll2004i iTremonti et, all 


2^p4 iSalim et al.l l200fl ~ Mannucci et, al.l 1201(1 120111 : 
Niinol I 2012 T ) . In this study, for our reference model, 


we adopt the galaxy stellar mass function and mass- 
dependent SFR proposed by iDrorv fe Alvared (l2008f) 
(hereafter DA08) as functions of redshift, and the 
redshift-dependent mass-metallicity relation from M08 is 
utilized. We also consider other models for comparison. 

3.1. Cosmic Star Formation Rate Density 

In DA08, the redshift evolution of the stellar mass 
function is presented o n the basis of the data spanning 
0 < 2 < 5 (IDrorv et al.l l2005h . It is assumed to have a 
Schechter form and is written as 


</>SMf(M*, 2 :) d-M* =(j) 0 [z) ^ a/DA 08 ( 


M* 


0 c), 

M* \ d M* 

6XP M^ A08 {z) ) Af 0 DA08 (^) ’ 

( 8 ) 

where the best-fitting parameterization is 

<j>o{z) = 0.0031 x (1 + z )~ 107 Mpc -3 dex -1 , (9a) 

fM» A08 ( z y 


log 


10 


V A 4 © 


= 11.35- 0.22 x ln(l + z), (9b) 


a DA08 = _ L3 _ 


(9c) 


The redshift dependence of SFR is also parameterized 
in DA08 as a function of stellar mass: 


M*(M,,*) = M“(aO 


M* 


exp - 


M DA 08 (,) 

M* \ 

M p A 08 (z) J ’ 


( 10 ) 
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with 


MP A 08 (z) = 2.7 x 10 10 x (1 + z ) 21 M q , (11a) 

/? DA08 = 0-6j ( llb ) 


while the fitting form of M®(z) is not given in DA08. In 
our study, we adopt the following functional form: 

M°(z) = M* N 0 15 ( 1 + *)“”“ exp ((3 m 5 z), (12) 

with M * N 0 15 = 1.183M 0 yr -1 , a N15 = 5.5 and / 3 N15 = 
—0.78, so as to fit the CSFRD, p*(z), of the data in 
DA08 shown in Figure [5] Note that CSFRD is obtained 
from 

/»00 

p*(z) = / dM*, (13) 

J o 

and the value of q 5 is determined to obtain p*( 0 ) = 
0.02M q yr -1 Mpc^ 3 , which is consistent with the 
values in recent studies of 0.017M Q yr -1 Mpc -3 
dHoriuchi et al.lf 2 (TTTh . O.O26M 0 yr -1 Mpc -3 (K13) and 
0.007-0.02M 0 yr -1 Mpc -3 (iMathew s et all 12 0141). 

So f ar, the CSFRD proposed bv iHonkins fo Beacoml 
(2006;) (hereafter HB06) has often been used for the cal¬ 
culation of the SRN spectrum. By performing a mod¬ 
ificatio n usi ng data of gamma ray bursts for z > 4, 
lYiiksel et al.l (1200811 obtained the parametric fit 


P* \Z) Po 


(l + z) c 


+ 


1 + 2 
C 


+ 


1 + 2 
B 

l/^HSce 


HB06 HB06 




(14) 


with a HB06 = 3.4, /3 HB06 = -0.3, y ™ 06 = -3.5, 
ryHB 06 _ —10 and po = O.O2M 0 yr -1 Mpc^ 3 . In 
this form, breaks are made at z\ = 1 and z 2 = 4 
by coefficients B = (1 + Zi ) 1 - qHB06 // 3HB06 and q = 
(1 + zx)^ )h (1 + z 2 ) l -P , respec¬ 

tively. While various observational results were compiled 
in HB06, there are uncertainties in the dust obscura¬ 
tion correction and the conversion from UV luminosity 
to SFR. In fact, most of the theoretical models for galaxy 
formation underpredict CSFRD at 2 > 1 and, according 
to K13, this inconsistency is due to these observational 
uncertainties. Therefore, for the SRN spectrum, we also 
examine CSFRD based on the semi-analytic model of 
galaxy formation presented in K13 as well as the obser¬ 
vational model in HB06. The CSFRD models considered 
in this paper are compared in Figure [ 2 ] We can see 
that our reference model of CSFRD based on DA08 lies 
between those in HB06 and K13 for z < 2, the region 
mainly responsible for SRNs. 

With the CSFRD models stated above, we obtain the 
total core-collapse rate, i?cc( 2 ), as 


d m 

IoTm® dM’ 


Rcc(z) = P*{z) X 


(15) 


adopting the Salpeter initial mass function (ip imf(-&0 c* 
M ~ 2 ' 35 ) with a mass range of O.1-1OOM 0 as in DA08. 



redshift 

Fig. 2.— CSFRD as a function of redshift. Dashed, solid and 
dotted lines correspond to the models in HB06, DA08 and K13, 
respectively. Plots are calculated from the data in Tables 1 and 2 
in DA08. 


Here M max and M m - m are the maximum and minimum 
masses of progenitors that end with a core collapse, re¬ 
spectively. Their values are given later in § [4] 

3.2. Metallicity Evolution 

The redshift dependence of the galaxy mass-metallicity 
relation was investigated in M08. The following analyti¬ 
cal form was adopted: 

12 + log 10 (O/H) = -0.0864 x (log 10 M* - log 10 M “ 08 ) 2 
+A^ 08 , (16) 

where the best fit values of log 10 M ^ 108 and K ^ 108 at dif¬ 
ferent redslrifts are listed in Table 5 of M08. While there 
are two choices for the values at 2 = 3.5 with differ¬ 
ent m odels of spectra l s ynthe sis, we adopt case a based 
on iBruzual fc Charlotl (12003f) . Using linear interpola¬ 
tion and extrapolation, we obtain the mass-metallicity 
relation for redshifts ranging from 0 to 5. We as¬ 
sume that the solar metallicity (Zq = 0.02) corresponds 
to the oxygen ab undance of 12 + log 10 (O/H) = 8.69 
ijAllende Prieto et al.il200ll ), i.e., 

logi ° = 12 + log 10 (O/H) - 8.69. (17) 


In our study, we assume that the stellar metallicity is 
identical to the value for its host galaxy for simplicity, 
whereas the metallicity vari ation inside a galaxy i s cur- 


rentlv being; discussed (e.£.. Levesaue et al. 

201 It iNiino 

2011 

Sanders et al. 2012: Taddia et al. 2013 

Niino et al. 

2014 

). We combine the mass-metallicity relation with the 


stellar mass function © and SFR (fTUl) of DA08 to de¬ 
rive the metallicity distribution function of progenitors 
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12+log 10 (0/H) 
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z/z 0 


12 + log 10 (0/H) 
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Fig. 3. — Normalized cumulative metallicity distribution function, which represents the fraction of progenitors with metallicity less than 
Z, for the models in DA08+M08 (left) and LN06 (right). The lines correspond, from bottom to top, to redshifts of z = 0, 1, 2, 3, 4 and 5. 


4>ZF(z,Z) as 


mass-metallicity relation was assumed to be 


V>zf (z,Z') d Z' 

rM,(z,Z) 

[0 


m,(ao)0smf(m:,*) dM: 


z) dM' 


(18) 


where M*(z, Z) is the stellar mass of a galaxy with metal¬ 
licity Z at redshift 2 calculated from Equations m and 
m The left-hand side of Equation (TTS1) , which is shown 
in Figure [3l represents the fraction of progenitors with 
metallicity less than Z. Hereafter, we refer to this model 
as DA08+M08. 

In this study, we also consider the model of me tallicity 
evolution proposed by lLanger fe Normanl (|2006f ) (here¬ 
after LN06). According to their model, the metallicity 
distribution function is given as 


r M,(z,Z) 


d Z’ = - 

Jo 


K<PsZf(M0 dM' 


K&iK) dM' 


T(a 


LN06 


2, M*(z, Z)/Mq N06 ) 


r(< 


V LN06 


2 ) 


(19) 


where T and T are the incomplete and com¬ 
plete gamma functions originating from the Scliechter 
form of their stellar mass function </>smf 6 (-^*) ^ 

(M* /Mq N06 ) “ LN06 exp(—M* /Mq N06 ), respectively. Note 
that, in Equation (fliil) . M*(j) smf is integrated instead of 
M t (f> smf in Equation m- However, M*</> smf is rea¬ 
sonable for our purpose because short-lived massive pro¬ 
genitors contribute to SRNs. Furthermore, in LN06, the 


M*(z, Z) 7 V-LNO 6 
M 0 LN06 



10 7 


LN06 ,qLN06 


( 20 ) 


where y LN06 
scaling as 


is related to the average cosmic metallicity 



= _ 7 LN06 


( 21 ) 


Taking the numbers in LN06, a LN06 = —1.16, /3 LN06 = 2, 
7 LNO 6 _ 0.15 and jF LN06 = 1 , we obtain the metallicity 
distribution function (Figure [3]). 


4. NEUTRINO EMISSION FROM STELLAR CORE 
COLLAPSE 

For the calculation of the SRN spectrum, we utilize the 
numerical results for the neutrino spectra emitted from 
various progenitors in the Supernova Neutrino Database 
(jNakazato et al.l 12013 ). The progenitor models in this 
data set a re co mp uted b y a Henyey-type stellar evolu¬ 
tion code (jlJmeda et a. 1 . 112012 ) and models with the ini¬ 
tial masses M = 13M 0 , 2OM 0 , 30M 0 and 50M 0 are 
provided. In our model, we set M m ; n = 10M 0 in Equa¬ 
tion (11511 for the minimum mass of progenitors that end 
with a core collapse for consistency with the evolution¬ 
ary calculation of the progenitors. We do not take into 
account the contribution due to electr o n-capture super¬ 
novae of low-mass progeni tors dLunardini fe Tamborral 
l2012t iMathews et al .1 [20 1 -111 . Note that recent theoret¬ 
ical s tudies suggest that their mass range is n arrow 
fe.g.. iPoelarends et al.l l2008t iDohertv et al.l l2015f) . On 
the other hand, the maximum mass of progenitors is 
M max = 1OOM 0 from the mass range of the initial 
mass function. Furthermore, in the Supernova Neu¬ 
trino Database, two different values of stellar metallicity 
(Z, = 0.02 and 0.004, i.e., the solar abundance and one- 
fifth of the solar abundance, respectively) are considered, 
and we consider the metallicity distribution for SRNs us¬ 
ing a threshold of Z = V0.02 x 0.004. Therefore, eight 
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Fig. 4. — Neutrino number spectra of supernova with 30 Mq, Z = 0.02 and shock revival times of Revive 
(solid) and 300 ms (dashed). The left, central and right panels correspond to v e , i'e and v x (= Vfx = = v T 


progenitor models are adopted in total. As shown in 
Equation m, their neutrino spectra, d N (M, Z,E' U )/ d E' v , 
are integrated in SRNs. 

4.1. Supernova Explosion and Shock Revival Time 

The physical factors that c ause the CCSN explosion 
are not well underst ood (e.g., IJankallSoTSD IKotake et al.1 
l2012tlBurrowsl[20l3ll . However, it is widely accepted that 
the shock wave launched by the bounce due to the nu¬ 
clear repulsion stalls once but is revived by some mech¬ 
anism. The process leading to shock revival is still a 
matter of debate. In the Supernova Neutrino Database, 
the shock revival time t rev ive is introduced as a parameter 
that reflects the unknown explosion mechanism. In this 
data set, under spherical symmetry, neutrino-radiation 
hydrodynamic simulations (lYam ada 1997: lYamada et all 
119991 : ISumivoshi et al.l 00051 ) and quas i-static evolu tion- 
ary calculations of neutrino diffusion (jSuzu kil FlOOl are 
used for the early and late phases of the supernova ex¬ 
plosion, respectively. Although the neutrino-radiation 
hydrodynamic simulations do not lead to a natural super¬ 
nova explosion, they are phenomenologically connected 
to quasi-static evolutionary calculations of neutrino dif¬ 
fusion,_based_on physical considerations (see Figure 14 
of I Na kazato e t all I2013D . The shock revival time af¬ 
ter the bounce has been estimated to be on the or der 
of 100 ms. For instance, according to iBelczvnski et al.l 
(|2012l) . it is preferably as short as 100-200 ms to ac¬ 
count for the observed mass distributions of neutron stars 
and b lack holes. On the other hand, lYamamoto et al.l 
(2013) suggested from a numerical study that the shock 
relaunch should be delayed until 300-400 ms to simul¬ 
taneously produce the appropriate explosion energy and 
nickel yields. Therefore, in our study, the shock wave 
is assumed to be revived at either t rev ive = 100, 200 or 
300 ms after the bounce. 

In Figure 01 we show the neutrino number spectra 
of supernovae with different shock revival times for the 
models with 30 Mq and Z = 0.02 in the Supernova Neu¬ 
trino Database. The total emission number and energy 
of supernova neutrinos depend on the shock revival time; 
they increase with Revive because more material accretes 
onto the collapsed core, releasing a huge amount of grav¬ 
itational potential energy. In particular, the difference is 
larger in the high-energy regime. This is because high- 
energy neutrinos are mainly emitted in the early phase 
when the proto-neutron star continues to be heated as a 


= 100 ms (dotted), 200 ms 
= u T ), respectively. 

result of mass accretion. After the shock revival, the ac¬ 
cretion stops and the mean energy of emitted neutrinos 
decreases, which is called the cooling phase. To estimate 
the uncertainty of the explosion mechanism, we calculate 
the spectrum of SRNs for models with different values of 
trevive (lNakazatoll201 3j). Note that the shock revival time 
may not be the same for all progenitors but be longer 
for more massive progenitors. We may underestimate 
the SRN flux using models with t rev i ve = 100 ms espe¬ 
cially for massive progenitors. Therefore, in this study, 
we regard the cases with f re vive = 300 and 100 ms as an 
indication of upper and lower bounds, respectively. 

4.2. Black Hole Formation and Nuclear Equation of 

State 

Some massive stars are thought to fail to pro¬ 
duce supernovae, leaving a bl a ck hole as a remnant 
(e.g.. iLi ebendor fer et al.ri2004 ISumivoshi et al.l 1200111 : 
IQ’Connor fc Otti I2011D . The failed supernovae also 
emit neutrinos from the bounce to black hole formation 
and, therefore, contribute to the overall flux of SRNs 
(|Lunardinil 120091 : iLien et all 120101 iNakazatol I2013D . In 
the Supernova Neutrino Database, the progenitor model 
with 30 Mq and Z = 0.004 is assumed to become a failed 
supernova because of its high core mass. Note that the 
models with Z = 0.004 have higher core mass than those 
with Z = 0.02 because the mass loss rate is larger for 
higher metallicity. On the other hand, the core mass is 
not monotonically related to the initial mass of progeni¬ 
tors because the mass loss rate is larger for a higher mass. 
Thus, the core mass of the model with the initial mass of 
30 Mq is the highest for the progenitors in the Supernova 
Neutrino Database. In contrast, progenitor models other 
than (Mj n it,Z) = (3OM 0 , 0.004) are assumed to be or¬ 
dinary supernovae. Combining this assumption with the 
Salpeter initial mass function ifiMpiM) and the metal¬ 
licity distribution function Z) of Equation (1181) or 

dm we obtain the fraction of black-hole-forming pro¬ 
genitors and calculate the spectrum of SRNs including 
failed supernovae. In Figure 01 the fraction of black-hole¬ 
forming progenitors is shown as a function of redshift for 
our models. At present, its val ue is hardly constrained by 
the observational data, while fHoriuchi et~aL (20141 ) esti- 
mat ed it to be 0.2-0.4 on the basis of another progenitor 
set (Wooslev et al.l [20021') . 

Although the neutrino signal from a failed supernov a 
is sensitive to the nuclear EOS (Sumivoshi et al.l [2006). 
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TABLE 1 

Numerical results for black hole formation of progenitor with (M,Z) = (30Mq, 0.004). 



£bh 

(E, e ) 

(Eo e ) 

(E, x ) 

Ei/ e ,tot 

Eu e ,tot 

Eu x ,tot 

Ev all, tot 

EOS 

(ms) 

(MeV) 

(MeV) 

(MeV) 

(10 52 erg) 

(10 52 erg) 

(10 52 erg) 

(10 53 erg) 

Shen 

842 

17.5 

21.7 

23.4 

9.49 

8.10 

4.00 

3.36 

LS(220 MeV) 

342 

12.5 

16.4 

22.3 

4.03 

2.87 

2.11 

1.53 


Note. — £bh is the time to black hole formation measured from the core bounce. The mean energy 
of the emitted i/» until black hole formation is denoted as (E Ui ) = E Ui , tot /E r u i , tot, where E Uii tot and 
Nv>i,tot are the total energy and number of neutrinos, respectively. u x stands for fi- and r-neutrinos 
and their anti-particles: E Ux = E Ufl = E&^ = E Ut = E„ T . ^i/ all ,tot is the total neutrino energy 
summed over all species. 


there is a result based on a single model by iShen et, ahl 
lfl998allbh in the Supernova Neutrino Database. To esti¬ 
mate the uncertainty of EOS, in this study, we compute 
the neutrino signal from the progenitor mode l with 30 Mq 
and Z = 0.004 using another EOS bv lLattimer &; Swestvl 
(j!991h (hereafter LS). While the LS EOS has three 
choices with different values of incompressibility K , we 
adopt the set with K = 220 MeV. The method of compu¬ 
tation is fully consistent with the original one in th e Su¬ 
pernova Neutrino Database (jNakazato et al.l[2013l) . In¬ 
cidentally, the EOS uncertainty of neutrino flux from 
ordinary C CSN has been evaluat ed for both the accre¬ 
tion ph ase (ISumivoshi et al.ll2005l ) and the cooling phase 
(|Snzukill2005l ) to be ~10%. which is smaller than the dis¬ 
persion due to the shock revival time in the Supernova 
Neutrino Database (~20-50%). 

We follow the neutrino emission until black hole for¬ 
mation with the neutrino-radiation hydrodynamic simu¬ 
lation. Note that the shock wave is not revived in this 
model and the mass accretion continues until black hole 
formation. Therefore, a larger amount of matter is ac¬ 
creted and more neutrinos are emitted compared with an 
ordinary supernova. Furthermore, the mean energy of 
neutrinos emitted from a failed supernova is higher than 
that from an ordinary supernova because the heating due 
to the accretion continues. It should also be noted that 



Fig. 5. — Fraction of black-hole-forming progenitors as a function 
of redshift. Dot-dashed and solid lines correspond to the models 
with the metallicity evolution of LN06 and DA08+M08, respec¬ 
tively. 


the neutrino emission from a failed supernova is rich in 
v e and i> e . This is because v e and v e are emitted more 
abundantly than v x (= = v T = v T ) from the ac¬ 

cretion of matter owing to the capture of electrons and 
positrons on nucleons. In contrast, during the cooling of 
a proto-neutron star, neutrinos of all species are emitted 
equivalently. 

In Figure [G] and Table [Q the results for the Shen EOS 
and LS EOS are compared. We find that the model with 
the Shen EOS produces more neutrinos than that with 
the LS EOS. This is for the following reason. The Shen 
EOS with the incompressibility of K = 281 MeV is stiffer 
than the LS EOS adopted in this study. Thus, the max¬ 
imum mass of neutron stars is higher and the amount 
of mass accretion is larger for the Shen EOS. In fact, 
the model with the Shen EOS takes a longer time to 
form a black hole (Table [lj. As a result, more poten¬ 
tial energy can be converted into the emission energy 
of neutrinos. These features are consistent with previ¬ 
ous studies (ISumivoshi et all 120061 iFischer et all 120091 : 
iNakazato et al.ll2010D . 


5. EXPECTED EVENT RATE AT SUPER-KAMIOKANDE 

In this section, we present the results for the spec¬ 
trum of SRNs and evaluate the event rate for Super- 
Kamiokande over 1 year. For this, we only consider the 
inverse /? decay reaction of electron antineutrinos, 

v e + p —> e + + n, (22) 

which is the most promising channel for detection. For 
Super-Kamiokande with a 22.5 kton fiducial volume, the 
number of target protons is set to N t = 1.5 x 10 33 . Note 
that, in the following, we do not take into account the 
detection efficiency for simplicity. Thus, using different 
normalization, our results are applicable to different de¬ 
tectors with the inverse (3 decay events. The event rate 
spectrum for positrons is written as 


dN e+ (E e+ ) 
d E e + 


N t a(E Pe ) 


d F(E Pe ) 
d E Pe 


(23) 


where cr(E Pe ) is the cross section for the inverse /3 de- 
ca y using the approx i matio n shown in Equation (25) 
of iStrumia fc Vissanil ( 2003 ). The positron energy is 
E e + = E Pe — Ac 2 with neutron-proton mass difference 
A. In the following, we choose a model with the CS- 
FRD of DA08, the metallicity distribution function of 
DA08+M08, a shock revival time of frevive = 200 ms and 
the Shen EOS as a reference model. 
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Fig. 6. — Neutrino number spectra for black hole formation with 30Mq, Z = 0.004 and Shen EOS (solid) and LS EOS (dotted). The 
left, central and right panels correspond to v (: , P e and v x (= = P fl — u T = u T ), respectively. 


5.1. Impact of Metallicity Evolution and Black Hole 
Formation 

We compare the results for the flux of SRNs, 
dF(Ep t .)/dEp t ., and the event rate spectrum, 
dN e +(E e +) / dE e +, for different metallicity evolu¬ 
tion models in Figure |TJ Here, we also show the results 
for fixed metallicity with Z = 0.02 and 0.004. The flux 
and event rate are higher for the case with Z = 0.004 
because th ey increase with the inclusion of failed 
super novae (|Lunardinil 120091 : ILien et all 1201 Ot INakazatol 
2013). In our model, the fraction of black-hole-forming 
progenitors for Z = 0.004 is fixed to ~0.14 while that 
for Z = 0.02 is 0. The accretion in failed supernovae 
continues until the mass of the proto-neutron star 
reaches the maximum mass of neutron stars. Therefore, 
more potential energy can be released from black hole 
formation than for ordinary CCSNe leaving neutron 
stars. Furthermore, the mean energy of neutrinos 
emitted from failed supernovae is higher than that for an 
ordinary CCSN because the neutrino-emitting surface 
of the core continues to be heated. Thus, the difference 
in the spectra for Z = 0.02 and 0.004 is larger in the 
high-energy regime. 

The enhancement of SRNs due to failed supernovae is 
larger for the normal mass hierarchy than for the inverted 
mass hierarchy because neutrinos from the accretion of 
matter dominate in the case of failed supernovae. As 
already stated, v e and v e are emitted more abundantly 
than v x (= = v T = v T ) from the accretion 

of matter. Nevertheless, for the inverted mass hierar¬ 
chy, all v e are converted to v^ or v T as in Equation 0. 
Meanwhile, for the normal mass hierarchy, the survival 
probability of v e is 0.68 as in Equation (0. On the other 
hand, the neutrino spectra of ordinary supernovae do not 
depend strongly on flavor because neutrinos of all species 
are emitted equivalently during the cooling of the proto¬ 
neutron star. Therefore, we can see the contribution of 
failed supernovae more clearly for the normal mass hier¬ 
archy. 

Models with the metallicity distribution functions of 
DA08+M08 and LN06 reside between the two fixed- 
metallicity cases with Z = 0.02 and 0.004 (Figure 0. 
The spectra of DA08+M08 and LN06 are similar. As 
shown in Tabled this is also the case for the event rates 
in various ranges of positron energy. In previous studies 
(lYtiksel h Kistlerll2012t INakazatol 20131) . the metallicity 
evolution of galaxies was not considered in the redshift 


dependence of the black hole formation rate. However, 
the black hole formation rate can be assessed using suit¬ 
able metallicity evolution models of galaxies. In fact, in 
our model, the uncertainty of the black hole formation 
rate in SRNs is small for the given mass and metallicity 
ranges of black-hole-forming progenitors. On the other 
hand, as already mentioned, the mass and metallicity de¬ 
pendences of the fate of progenitors are still unclear. For 
the normal mass hierarchy, the event rate in the positron 
energy range of 10-26 MeV for the model with a fixed 
Z of 0.004 is ~40% larger than that when Z = 0.02. 
Thus, when the fract ion of black -hole-forming progeni¬ 
tors is ~0.3 (jHoriuchi et a,l.lf20~i~il ). the event rate would 
become ^80% larger than that for the model with Z fixed 
to 0.02. In contrast, for the inverted mass hierarchy, the 
enhancement of the event rate is ~15% for Z = 0.004 
and, therefore, it would be ~30% when the fraction of 
black-hole-forming progenitors is ~0.3. 

Now we move on to the EOS dependence of black hole 
formation while the above investigations in this subsec¬ 
tion were based on the models with the Shen EOS. In 
Figure [HJ we show the results obtained with the Shen 
EOS and LS EOS for fixed metallicity with Z = 0.004 
for the purpose of co mparison. We can see that, as an¬ 
ticipated from § 14.21 the flux and event rate are lower 
for the LS EOS. In general, a stiffer EOS can result 
in higher SR N flu x, which is consistent with a previous 
study (iLunardinilI2009D . In particular, for the LS EOS, 
the flux of the model with Z = 0.004 does not signifi¬ 
cantly differ from that with Z = 0.02, which does not 
include failed supernovae. The contribution of failed su¬ 
pernovae is also not clear for the event rates in various 
ranges of positron energy (Table 0. In conclusion, the 
enhancement of SRNs due to failed supernovae is signif¬ 
icant for the case with the normal mass hierarchy and a 
stiff EOS. 

5.2. Dependence of Cosmic Star Formation Rate 
Density 

The results for the SRN flux and event rate spectrum 
for various CSFRD models are shown in Figure 0 A 
model with a higher CSFRD has more progenitors of 
both ordinary and failed supernovae as in Equation (1151) . 
We can recognize that the difference is clearer for the 
low-energy regime, in contrast to the models in Figure 0 
This is because, reflecting the fact that the uncertainty 
in CSRFD is larger for a higher redshift, the difference 
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Normal Inverted 



E e + [MeV] E e +[MeV] 


Fig. 7. — Fluxes of SRNs (upper panels) and event rate spectra in Super-Kamiokande over 1 year (lower panels ) obtained using models 
with CSFRD of DA08, shock revival time of t rev ive — 200 ms and Shen EOS. The left and right panels show the results for the normal and 
inverted mass hierarchies, respectively. Solid and dot-dashed lines correspond to models with the metallicity evolutions of DA08+M08 and 
LN06, respectively, while other lines denote the results for fixed metallicity with Z = 0.02 (dotted) and 0.004 (dashed). 


Normal Inverted 




E e + [MeV] 


E e +[MeV] 


Fig. 8. — Same as Figure [7] but for different EOSs. Dashed and dotted lines correspond to models with the Shen EOS and LS EOS, 
respectively, with the metallicity fixed to Z — 0.004. The contribution of black-hole-forming failed supernovae is extracted for the Shen 
EOS (dot-short-dashed) and LS EOS (dot-long-dashed). Solid lines denote the results for fixed metallicity with Z = 0.02, where failed 
supernovae are not included. 
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TABLE 2 

SRN EVENT RATES IN VARIOUS RANGES OF POSITRON ENERGY IN SUPER-KAMIOKANDE OVER 1 YEAR 
(I.E., PER 22.5 KTON YEAR) FOR MODELS WITH CSFRD OF DA08 AND SHOCK REVIVAL TIME OF 

irevive = 200 MS. 


metallicity evolution 

EOS for BH 

Normal mass 

hierarchy 

Inverted mass 

hierarchy 

18-26 

10-18 

10-26 IVIeV 

18-26 

10-18 

10-26 IVIeV 

Z = 0.02 


0.227 

0.549 

0.776 

0.301 

0.640 

0.941 

DA08+M08 

Shen 

0.274 

0.604 

0.879 

0.326 

0.660 

0.986 

LN06 

Shen 

0.288 

0.625 

0.912 

0.334 

0.669 

1.003 

Z = 0.004 

Shen 

0.387 

0.714 

1.100 

0.378 

0.691 

1.069 

Z = 0.02 


0.227 

0.549 

076 

0.301 

0.640 

0.941 

DA08+M08 

LS 

0.233 

0.554 

0.787 

0.308 

0.640 

0.948 

LN06 

LS 

0.235 

0.557 

0.791 

0.311 

0.642 

0.953 

Z = 0.004 

LS 

0.247 

0.563 

0.810 

0.322 

0.632 

0.954 


Normal 


Inverted 




E e + [MeV] E e +[MeV] 

Fig. 9. — Same as Figure[7]but for the models with metallicity evolution of DA08+M08, shock revival time of t rev ive = 200 ms and Shen 
EOS. Dashed, solid and dotted lines correspond to models with the CSFRD of HB06, DA08 and K13, respectively. 


Normal Inverted 




Fig. 10.— Total fluxes of SRNs (solid) and contributions from various redshift ranges for the reference model. The lines except for the 
solid line correspond, from top to bottom, to the redshift ranges 0<z<l, l<z<2, 2<z<3, 3<z<4 and 4 < z < 5, for 
E u > 10 MeV. The left and right panels show the cases for normal and inverted mass hierarchies, respectively. 
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Normal 



E e *[MeV] 


Inverted 



E e *[MeV] 


Fig. 11.— Same as Figure[7]but for models with CSFRD of DA08 and metallicity evolution of DA08+M08 and Shen EOS. Dotted, solid 
and dashed lines correspond to models with shock revival times of t rev ive = 100, 200 and 300 ms, respectively. 


in CSFRD among the models adopted here is large for 
the redshifts with z > 0.5. Neutrinos from the dis¬ 
tant Universe accumulate in the low-energy region of the 
SRN spectrum because they are redshifted by the cos¬ 
mic expansion and their energy is reduced by a factor of 
(1 + 2 ) . For the SRN flux of the reference model, con¬ 

tributions from various redshift ranges are indicated in 
Figure [TUI We can see that, a s discussed in previou s 
works (e.g., lAndo fe Satol 12004 IMathews et all 120141 . 
the considerable flux of low-energy SRNs is attributed 
to the range 1 < 2 < 2 while the dominant flux origi¬ 
nates from the local Universe (0 < z < 1) for neutrino 
energy E v > 10 MeV. 

5.3. Dependence of Shock Revival Time 

In Figure fill we show the spectra obtained from the 
models with different values of the shock revival time, 
which is introduced as a para meter reflecting the st ill un¬ 
known explosion mechanism dNakazato et al.l[2013lf . The 
models with a longer shock revival time have a higher 
event rate because, as already stated in § 14.11 the total 
emission number and energy of supernova neutrinos in¬ 
crease with the shock revival time. The increase in flux 
is clear for the high-energy regime; this trend is simi¬ 
lar to that observed with the inclusion of failed super- 
no vae ( § 15.11 and Figure 0. Nevertheless, as reported 
by iNakazatol (l2013ll . the hardening of the SRN spec¬ 
trum is not significant because a considerable fraction 
of neutrinos is emitted after the shock revival, where a 
proto-neutron star is not heated and the mean energy of 
neutrinos gradually decreases. Note that, the increase in 
flux (Figure 1111) is larger for the normal mass hierarchy 
because u e and v e are abundant in the accretion phase, 
which is longer in the case of later shock revival. 


5.4. Summary of Uncertainties 

Here, we summarize the uncertainties in the SRN spec¬ 
trum due to the EOS for black hole formation, CSFRD 
and shock revival time t lev - lve , adopting the models with 
DA08+M08 for metallicity evolution, whose variance was 
shown to be small in §0 In Table O the event rates in 
various ranges of positron energy are listed for all com¬ 
binations. We can see that the uncertainty of CSFRD 
is relatively larger than those of the EOS and f rev ive for 
events in a low-energy range (10 MeV < E e + < 18 MeV), 
whereas they are comparable for events in a high-energy 
range (18 MeV < E e + < 26 MeV). This is because the 
SRN flux dependences on the EOS and f rev ive are mi¬ 
nor for low energies. This tendency is especially clear in 
the case of the inverted mass hierarchy, where the de¬ 
pendences on the EOS and i rev ive are again minor. In 
Figure m we show the spectra of the reference model 
and models with the maximum and minimum values of 
the SRN event rate given in Table [3] We can recognize 
that the uncertainty in the SRN flux is a factor of two 
or three across the full energy range shown in Figure [TUI 
As already described, the difference in the low-energy 
spectrum mainly originates from CSFRD (see Figure El), 
whereas the choice of EOS and t rev ive also affects the 
high-energy spectrum (see Figures l8l and fTTl) . 

We compare the results for the normal and inverted 
mass hierarchies in Figure [13] The difference is clear for 
model with minimum values of SRN event rate, where a 
short shock revival time (frevive = 100 ms) and LS EOS 
are adopted. In models with short shock revival time 
(see Figure [4]) and LS EOS (see Figure [6]), the difference 
in the spectra for v e and v x before the neutrino oscilla¬ 
tion is larger in the high-energy regime. This is for the 
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TABLE 3 

SRN EVENT RATES IN VARIOUS RANGES OF POSITRON ENERGY IN SUPER-IvAMIOKANDE OVER 1 YEAR (i.E., PER 
22.5 KTON YEAR) FOR MODELS WITH METALLICITY EVOLUTION OF DA08+M08. 


CSFRD 

^revive 

EOS for BH 

Normal mass 

hierarchy 

Inverted mass 

hierarchy 

Figure 1121 

18-26 

10-18 

10-26 IVLeV 

18-26 

10-18 

1U-26 IVieV 

HB06 

100 ms 

Shen 

0.286 

0.704 

0.990 

0.375 

0.832 

1.207 




LS 

0.227 

0.635 

0.863 

0.351 

0.806 

1.156 



200 ms 

Shen 

0.361 

0.833 

1.193 

0.429 

0.920 

1.349 




LS 

0.302 

0.764 

1.066 

0.404 

0.893 

1.297 



300 ms 

Shen 

0.432 

0.938 

1.370 

0.463 

0.967 

1.431 

Maximum 



LS 

0.374 

0.869 

1.242 

0.439 

0.941 

1.379 


DA08 

100 ms 

Shen 

0.219 

0.515 

0.734 

0.286 

0.598 

0.885 




LS 

0.178 

0.464 

0.642 

0.269 

0.578 

0.847 



200 ms 

Shen 

0.274 

0.604 

0.879 

0.326 

0.660 

0.986 

Reference 



LS 

0.233 

0.554 

0.787 

0.308 

0.640 

0.948 



300 ms 

Shen 

0.326 

0.677 

1.003 

0.350 

0.694 

1.044 




LS 

0.285 

0.627 

0.911 

0.333 

0.674 

1.007 


K13 

100 ms 

Shen 

0.203 

0.443 

0.645 

0.264 

0.505 

0.769 




LS 

0.171 

0.410 

0.581 

0.252 

0.492 

0.744 

Minimum 


200 ms 

Shen 

0.252 

0.514 

0.767 

0.298 

0.554 

0.853 




LS 

0.221 

0.482 

0.703 

0.286 

0.542 

0.827 



300 ms 

Shen 

0.298 

0.570 

0.868 

0.319 

0.580 

0.899 




LS 

0.266 

0.537 

0.804 

0.306 

0.568 

0.874 



Normal Inverted 





Fig. 12.— Same as Figure[7]but for reference model (solid) and models with maximum (dashed) and minimum (dotted) values of SRN 
event rate among models with metallicity evolution of DA08+M08. See Table [3] for the parameter sets of the maximum and minimum 
models. 


following reason. The emission of high-energy v e and v e 
is suppressed owing to dense accreting matter in early 
stage of the accretion phase (<100 ms after the bounce). 
The average ener gy of v x is higher eve n in the early stage 
(see Figure 14 of lNakazato et al.1 120131 ) because a mean 
free path of v x is longer than those of v e and v e . On the 
other hand, high-energy v e and v e are more efficiently 
emitted than v x later on the accretion phase (>200 ms 
after the bounce) and the difference in the total emission 
number of v e and v x gets smaller. Thus the difference 


in the SRN spectra for the normal and inverted mass 
hierarchies is clear for models with a short mass accre¬ 
tion. Note that LS EOS model has a shorter time to 
black hole formation, which corresponds to the duration 
of mass accretion, than Shen EOS model. 

6. CONCLUSION AND DISCUSSION 

In this paper, we have studied the supernova relic neu¬ 
trino (SRN) spectrum and event rate involving black- 
hole-forming failed supernovae. To evaluate the contri¬ 
bution of failed supernovae to SRNs, cosmic metallicity 
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Fig. 13.— Event rate spectra in Super-Kamiokande over 1 year 
obtained using reference model (solid) and models with maximum 
(dashed) and minimum (dotted) values of SRN event rate among 
models with metallicity evolution of DA08+M08. Thick and thin 
lines correspond to the results for the normal and inverted mass 
hierarchies, respectively. 


evolution is important as well as the mass and metallic¬ 
ity dependences of the fate of progenitors. In this study, 
impact of metallicity evolution on SRNs has been investi¬ 
gated for the first time. We have determined the redshift 
evolution of the metallicity distribution function of pro¬ 
genitors from observed galaxy stellar mass function and 
empirical relation of galaxy mass and metallicity. As a 
result, the uncertainty of the black hole formation rate in 
our model spectra of SRNs is small for the given mass and 
metallicity ranges of black-hole-forming progenitors. Our 
reference model of SRNs has a lower flux compared with 
previous studies, where a high er cosm i c sta r formation 
rate density (CSFRD) (jHonkins fe Beacomll2006l) and a 
higher mean e nergy of emitted neutrinos were assumed 
(e.g., lBeacomll201(ll ). 

We have also investigated the dependences of SRNs 
on the CSFRD, shock revival time t rev ive and equation 
of state (EOS). The shock revival time is introduced as 
a parameter that should depend on the still unknown 
explosion mechanism of core collapse supernovae. The 
EOS dependence has been considered for failed super¬ 
novae, whose collapse dynamics and neutrino emission 
are certainly affected, using the new result on a failed 
supernova with a different EOS. It has been found that 
the differences with regard to t rev ive and the EOS are 
clear for the high-energy regime, especially in the case of 
the normal mass hierarchy. In contrast, the low-energy 
spectrum of SRNs is mainly determined by CSFRD for 
both mass hierarchies. Therefore, the detection of low- 
energy neutrinos is mandatory so as to probe the cosmic 
star formation history by SRN observation. 

Although numerous noise events prevent the observa¬ 
tion at low energies, considerable effort is now being de¬ 
voted to reducing them. Photomultipliers are used to 
observe Cherenkov photons (in a water Cherenkov type 


detector such as Super-Kamiokande) or scintillation pho¬ 
tons (in a liquid-scintillator-type detector such as Kam- 
LAND) emitted from positrons produced in the inverse 
beta decay (l22l) . Meanwhile, the neutrons produced in 
Equation (l2fil) are captured on protons and emit de¬ 
layed gamma rays (2.2 MeV). Therefore, v e events are 
identified by detecting 2.2 MeV gamma rays in Kam- 
LAND and recently in Super-Kamiokande dZhang et al.1 
I2015D . Distinguishing D e events from other background 
events without neutron emission such as invisible muon 
events is called neutron tagging. Furthermore, neu¬ 
tron tagging will be more efficient with gadolinium- 
loaded water because of the high neutron capture rate 
on gadoli nium and the high-energy (8 MeV) gamma-ray 
emission (Be acom fc Vagini 12004) . Actually, gadolin¬ 
ium loading is planned in Super-Kamiokande. This will 
allow the detection of low-energy events in the range 
of 10 MeV < E e + < 18 MeV investigated in our 
study as a regime where the event number is mainly 
determined by the CSFRD. In fact, assuming 90% o f 
the neutron-tagging efficiency dWatanabe et al.l I2009D , 
Super-Kamiokande will detect 4-9 low-energy events over 
10 years depending on CSFRD. We hope that Super- 
Kamiokande with gadolinium-loaded water will become 
a powerful detector of SRNs, opening a new door to ex¬ 
ploring the history of the Universe. 

The authors are grateful to Makoto Sakuda and 
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Physics (RCNP) in Osaka University. This work was par¬ 
tially supported by Grants-in-Aids for the Scientific Re¬ 
search (No. 24105008, No. 26104006, No. 26105515 and 
No. 26870615) from MEXT in Japan. 



Fig. 14.- Fluxes of SRNs obtained using reference model for the 
case without neutrino oscillation. Dashed, solid and dotted lines 
correspond to i/ e , P e and v x (= u fl = = v T = P r ), respectively. 


APPENDIX 

SRN FLUX DATA TABLE 

For general use in any research for astronomy, astrophysics, and physics, the numerical data of SRN flux in this 
study are publicly available on the Web at 

http://asphwww.ph.noda.tus.ac.jp/srn/ 

We provide not only the reference model but also models with maximum and minimum values of SRN event rate 
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among models with metallicity evolution of DA08+M08. Moreover, for a broad application, fluxes of SRNs with i/ e , 
P e , and v x for the case without neutrino oscillation, which are shown in Figure IT4l are also available. 
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